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Abstract 

We apply the relative weights method (arXiv: 1209.5697) to determine the effective Polyakov line action 
for SU(2) lattice gauge theory in the confined phase, at lattice coupling j3 = 2.2 and N t =4 lattice spacings in 
the time direction. The effective action turns out to be bilinear in the fundamental representation Polyakov 
line variables, with a rather simple expression for the finite range kernel. The validity of this action is tested 
by computing Polyakov line correlators, via Monte Carlo simulation, in both the effective action and the 
underlying lattice theory. It is found that the correlators in each theory are in very close agreement. 
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I. INTRODUCTION 



The Polyakov line action (PLA) is an action obtained from lattice gauge theory when all degrees 
of freedom are integrated out, under the constraint that the Polyakov line holonomies are held 
fixed. There are some indications [1-4] that the sign problem in this theory, at non-zero chemical 
potential, may be more tractable than the sign problem in the underlying lattice gauge theory (for 
a review, cf. [5]), and if so it could provide us with a new tool for investigating the QCD phase 
diagram. It is fairly straightforward, given the PLA at chemical potential /i = 0, to introduce a 
non-zero chemical potential, as we discuss in section V. The problem we address here is how to 
extract the PLA from the underlying lattice gauge theory at /i = 0. 

This article is a follow-up to ref. [6], which presented a novel "relative weights" technique 
for deriving the PLA, based on a method used previously in studies of the Yang-Mills vacuum 
wavefunctional [7]. The method was tested at strong couplings, where the answer is known, and a 
conjecture for the action at a weaker coupling, for SU(2) pure gauge theory at /3 = 2.2 and inverse 
temperature N t =4 lattice spacings, was presented. This conjecture was based, however, on a study 
limited to fairly atypical regions of field configuration space. Below we will apply the method in 
the region expected to dominate the path integral, and a rather different (and in fact simpler) action 
from the one conjectured in ref. [6] emerges. As a crucial test of the derived PLA, we compute 
the two-point Polyakov line correlator from Monte Carlo simulations of both the PLA and the 
underlying gauge theory. These correlators will be seen to agree quite accurately. Our effective 
PLA turns out to be only bilinear in the Polyakov line variables, with a simple expression for the 
finite range kernel. 

There have been a number of previous attempts to derive the PLA from lattice gauge theory at 
finite temperature. These include strong-coupling expansions [2], the Inverse Monte Carlo method 
[8], and the demon approach [9]. All of these methods generate effective Polyakov line actions of 
varying degrees of complexity. We believe, however, that an accurate agreement of Polyakov line 
correlators in the confined phase, computed in the effective and underlying lattice gauge theories, 
has not been demonstrated in any of the previous studies, at least not beyond two or three lattice 
spacings in Polyakov line separation. 

It should also be mentioned that there are a number of studies which are concerned with de- 
ducing the Polyakov line potential, with particular application to the deconfinement transition, cf. 
[10] and references therein. There have also been efforts, e.g. [11], to express the fermion deter- 
minant in terms of a potential involving Polyakov lines. These studies do not arrive at a full PLA 
as defined above, and hence their focus is somewhat different from ours. 

Our article is organized as follows: An improved version of the relative weights method is 
presented in section II below. The technique is applied to pure SU(2) gauge theory in section III, 
again at /3 = 2.2 and N t = 4, and the Polyakov line correlators of the derived PLA are compared 
to those of lattice gauge theory. Application to a gauge-Higgs theory, with a scalar matter field 
explicitly breaking global Z2 center symmetry, is presented in section IV. Section V contains our 
conclusions. The extension of our method to gauge theories with dynamical fermions is discussed 
in an appendix. 

H. THE RELATIVE WEIGHTS METHOD 

The relative weights method, as applied to deriving the PLA Sp, was introduced in ref. [6]. 
The technique is particularly well adapted to computing path (or "directional") derivatives of the 
effective action Sp in the space of all Polyakov line configurations. A single configuration {L^} 
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is a point in this space, where we specify the group-valued Polyakov line holonomies U x at each 
spatial point jc in a three-dimensional volume. Let {U X (X)} be a path through this space of config- 
urations, where A parametrizes the path. The relative weights method computes the path derivative 
dSp[U x (X)]/dX at some given A = Xq. In this section we present a variant of the relative weights 
approach which, while equivalent to the original method of [6], is numerically more efficient. 



In order to minimize minus signs later on, we adopt the convention that the Boltzmann weight 
is proportional to exp [+£/>]. Let Sl be the lattice gauge action on an L 3 x iV f volume with coupling 
/3 for the Wilson action. Sl may contain pseudofermion or bosonic matter degrees of freedom, 
collectively denoted by 0. It is convenient to go to temporal gauge, so that all timelike link 
variables are set to the unit matrix except on a timeslice at t = 0. Then the PLA Sp is defined as 



exp 



Sp[U x ]}= [ DU (x,0)DU k D<j> \l\8[U x -U (x,0)}\e SL 



(1) 



Because of the residual Uq(x,0) — > g(x)Uo(x,0)g^ (x) symmetry in temporal gauge, it follows that 
Sp can only depend on the eigenvalues of the U x matrices. 



While the functional integration in (1) can only be carried out in special cases, e.g. via strong 
coupling and hopping parameter expansions valid a certain range of parameters, the ratio (or "rel- 
ative weights") exp[Sp [£/£]]/ exp[5p[f/^']] evaluated at nearby configurations U' X1 U' x ' is calculable 
numerically. This fact enables us to compute path derivatives of Sp. 



Let us consider a set of M Polyakov line configurations 

{{C/j n) ,allx}, ti=1,2,...,m} 
corresponding to values of the path parameter 

K = Xq + ( n - 1 ) AA , n = 1,2,...,M 



and define 



S L m) [U, <j>] = S L U (x, 0) = UP ,U k {x,t),${x,t) 



(2) 



(3) 



(4) 



to be the lattice action in temporal gauge with the timelike links at t = fixed to the m-th member 
of the set (2). We also define 



AS£" +1) = S P [u(' n+ ^] -S P [U^] 



DU k D(j) e" 



Jm) 



(5) 



From (1), we have 



exp[As£" +1) 



exp 




exp 





(m+1) 



fPU k D$ e S - 
fDU k D(j) e s i' 
jDU k D(j) exp 



Jm) 



fDU k D(j) e s L 



exp 



AS 



m+l) 



(6) 



where (...) m indicates that the expectation value is taken from ensembles with Boltzmann factor 
exp [S^]/Z m . For sufficiently small AX, 

ds P [u x (x)} as { ; i+1) 



dX 



AX 

AT log( 



AS 



(m+l) 



(7) 



and these should closely agree, for all m < M, with the derivative evaluated at the central value 
of X = Xq. We can then improve our estimate by making use of all M configurations, taking the 
average of derivatives 



1 1 



f dS P [U x {X)] \ 
V dX ) x=h 'AXM 



M-l 



I E lo g(( ex p 



m=l 



AS 



(m+l) 



(8) 



The question then becomes which point {U x (Xq}} in configuration space should be chosen for 
the computation, and which directional derivatives dSp/dX at this point should computed, in order 
to deduce Sp. It is possible that the choice is not very important, and that Sp is well approximated 
by the same simple expression everywhere in configuration space. However, if this is not the case, 
then calculating path derivatives in some very atypical corner of configuration space may lead to 
an approximate answer for Sp which may be correct in that particular corner, but misleading in the 
bulk of configuration space. 

In ref. [6] the path derivatives were computed using three types of sets (2) for SU(2) lattice 
gauge theory. These were (i) Polyakov lines which were constant in space, and the X parame- 
ter was the amplitude of P x = jTt[U x ] = X; (ii) Polyakov lines which consisted of small plane 
wave fluctuations around a constant background, P x = Pq + X cos(k ■ x) with X <C Po; and (iii) 
Polyakov lines in which P x varied as P x = X cos(k ■ x). Compared to thermalized timeslice con- 
figurations £/o(*, 0) generated in a normal lattice Monte Carlo simulation, such configurations are 
very atypical. In a thermalized configuration in the confined phase, the Fourier components of the 
configuration are all of 0(1 /\/Vi), where V3 = L 3 is the lattice volume of the D=3 dimensional 
timeslice, whereas in the special configurations just mentioned, one Fourier component (which 
may be the zero mode) is of 0(1). For computing the PLA at a strong lattice coupling, where 
this action can be evaluated via a strong coupling expansion, the atypical nature of the constant + 
plane wave, or pure plane wave configurations did not seem important, and the PLA deduced from 
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the relative weights data was a close match to the known result. There was no similar result to 
compare to at /3 = 2.2, A// = 4, and although an expression for Sp matching the results for dSp/dX 
was deduced from fitting the data, there is a concern that this expression might only be valid in the 
special region of configuration space where it was derived. 

We will investigate the action in more typical regions of configuration space in the next section. 



III. DERIVATIVES OF S P IN A THERMALIZED BACKGROUND 

The goal is to use the expression for the path derivative (8) to determine Sp. However, if Sp 
does not have a simple form everywhere in configuration space (and it may not), then it is at least 
required that we have a fairly accurate approximation to Sp in the region which is important for the 
computation of observables, i.e. the region occupied by typical fhermalized configurations {U x }. 
A set of timelike link configurations {Uq(x,0)} on the t = timeslice, generated by a numerical 
simulation of the underlying lattice gauge theory, is a sample of such configurations. Let us define, 
for the SU(2) gauge group that we will consider here, 



where the sum runs over all wavevectors k on a cubic lattice of volume L , and = a = —b k 
are real-valued. Then we may consider calculating numerically, by the relative weights method, 
derivatives with respect to the Fourier components 



in a background in which all other Fourier components are drawn from a thermalized configura- 
tion. By calculating the derivative for some range of a, it is possible to extrapolate to small a 
of order 1 / \/Vi, which is the typical magnitude of a Fourier component in thermalized config- 
urations. It is sufficient in practice to concentrate on the coefficients of the cosine terms in the 
Fourier expansion, since the sine terms give similar results. We then try to reconstruct Sp from 
this information. 

There are potentially two obstacles to this approach. First, there are as many independent 
Fourier components as there are lattice sites in the V3 volume, and this is too many to calculate 
in practice. Secondly, it might be that the results are strongly dependent on the particular ther- 
malized background which is used. Concerning the first obstacle, this will not be a problem if the 
derivatives with respect to have a simple dependence on the lattice momentum 



which can be deduced from a small sample of all possible components. As for the second obstacle, 
this is not a problem if it turns out that the dependence of the final results on the particular choice 
of thermalized configuration is very weak. 




(9) 




(10) 




(11) 
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A. Deriving the effective PLA 



The first step in the extrapolation to small a is to run a standard lattice Monte Carlo, stop at 
some thermalized configuration, and calculate all the Polyakov line holonomies, 

& x = Uo(x,l)U {x,2)..-Uo{x,N t ) 

= d4(ni,n 2 ,n 3 )l + id(ni,ri2,n3) ■ a , (12) 

where x = («i, n 2 , "3)> m that configuration. Define W(x) = d^{x). We pick a particular wavevector 
k which is specified by three integers (mi , m 2 , m 3 ) with corresponding wavenumber components 

h = -j-m , (13) 

and set the coefficient of cos(fc-jc), in the sine-cosine expansion of W(x), to zero. Denote the 
modified array, with the cos(fc-jc) term removed, as W'(x). 

Next, construct a set of M = 20 configurations with: 

= fl W cos(kmi + k 2 n 2 + k 3 n 3 ) + (1 - a - 8)W'(x) 

a W = a+(n--(M+l)^y/L 3 , n=l,2,..,M, (14) 

where 7 = L 3 Aa is a constant chosen to be as small as possible, but still large enough to get some 
spread in the data. Typically 7 ps 0.5. 

The factor 1 — a — 8 in (14) is introduced in order to keep P x , with rare exceptions, inside 
the range [—1 : 1]. Ideally one would like to leave all {a^,b^} in the thermalized configuration 

unaltered, apart from the mode with k' = k, i.e. P x = akCosk-x + W^x). At finite a, however, 
this has the disadvantage that at many sites \P X \ > 1. To see this, note that W(x), from which 
W'(x) is derived, may come close to the limits ±1 at some sites, and at these sites the additional 
contribution a/ c cos(k ■ x) may put the sum outside the allowed range by as much as a. Moreover, 
by removing the cos(fc-Jt) mode, W'(x) may already lie outside the range [—1:1] at some sites; 

(n) 

this is especially true for k = 0. We must also allow for the fact that half of the {a k } are slightly 
greater than a. For this reason we reduce the amplitude of the added thermalized configuration 
by a factor of 1 — a — 8 (in our simulations we found 8 = .04 sufficient). In the exceptional cases 

(n) 

where P x still lies outside the allowed range, it is truncated to the nearest limit, i.e. ±1. 

in) 

We then construct the SU(2) variables, at each site, which have 2P X ' as the trace 

t/i n) =pi n) l + ^(jc)-<7, (15) 

where, to insure unitarity, 



/ 1 ~ ( P *" ) 2 (16) 
The calculation of {dSp/ da^) a proceeds as described above. For the choice of U^ n > given above, 
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L=24 data at a=0.05 




FIG. 1. Derivatives of the PLA L 3 dSp/da/ c evaluated at = a = 0.05, vs. lattice momenta k^. Also 
shown is a linear best fit to the data at &l > 0.7. 



L=24 data with a rescaling 
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FIG. 2. Derivatives L 3 (dSp/dai c ) a divided by a, vs. lattice momenta k L , for a = 0.05,0.10,0.15,0.20. It 
is clear that the derivatives of Sp depend linearly on a. 



it is easy to see that for a lattice of extension L in the spatial directions 

The results for the derivatives are found to depend only weakly on the choice of thermalized time 
slice (12) generated in an ordinary Monte Carlo run. The dependence is most pronounced at small 

with the variance on the order of 2%. In practice we have averaged our results for dSp/da^ 
over eighty independent time slices. 

In Fig. 1 we display our results for the L~^dSp/da^ vs. k L at a = 0.05, and lattice spatial 
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extension L = 24. (Note that, apart from Figs. 5 and 7, all our figures show data derived at an 
L = 24 extension.) The underlying SU(2) lattice gauge theory is defined as a Wilson action on a 
periodic 24 3 x 4 volume, at the coupling /3 = 2.2. The calculations were made, in this case, at 
lattice momenta with components ki = Inrm/L, with the following {m\m2m-s) triplets: 



(000), (100), (110), (111), (200), (210), (211), (300), (311), (320), 
(400) , (322) , (42 1 ) , (430) , (333) , (433) , (443) , (444) , (554) , (654) , 
(655), (665), (766), (777), (887), (988), (998), (1099), (10 10 10)) . (18) 



On this plot the data point displayed at kt = is a factor of two smaller than the actual data value; 
this was done for reasons to be explained shortly. 

The striking thing about this data is that, for lattice momenta kt > 0.7, the data points clearly 
fall on a straight line. The second fact is that the data is linearly proportional to a at these small 
a values, as we see in Fig. 2. In this figure we divide dSp/da^ by a, at a = 0.05,0.10,0.15,0.20, 
and find that the data points coincide. The linearity of the derivative w.r.t. a k implies that the action 
itself is quadratic in these variables, leading to a simple bilinear form 



Sp = ^c l y £P x l -2c 2 Y,PxQ(x-y)Py , 

1 x xy 

where 

Q(x-y) = ^Q(kL)e' k - {x - y) . (20) 



4 M*-y) 



This leads to derivatives 



J_ ( dS P [U x (a k )] \ = \ tt( 5 Cl " 2c 20(^)) k L ± 
L 3 V da k J ak=a | 2a(ici-2c 2 G(0)) k L = 



(21) 



2' 

The relative factor of two in the kt = and ki > cases is due to the fact that 1 = L 3 , while 
£ x cos 2 (&- x) = \L 3 . The kL > data should extrapolate, as kt — > 0, to a value which is half the 
result at kt = 0, which is why we have divided the derivative at kt = by a factor of 2, when 
displaying these values on Figs. 1 and 2. The constants c\ and c% are obtained from a linear fit to 
the data at kt > 0.7, as shown in Fig. I. 1 

It is clear that for kt > 0.7, the fc-space kernel is Q(k£) = kt- If this were true at all kt, then we 



would have Q = y — in position space, where is the lattice Laplacian. However, a kernel 
of this kind is infinite range, which would violate one of the assumptions of the Svetitisky-Yaffe 
analysis (cf. [12]). In any case, Q(k£) deviates from linearity at small momentum. We therefore 



1 In practice we fit the data for each a, at k L > 0.7, to the form A(a) —B(a)k L . We then fit A(a), B(a) to straight 
lines, and the constants c\, C2 are extracted from the slopes, i.e. dA/da = jCi, and dB/da = 2c2- The choice of 
0.7 as the lower limit is a potential source of systematic error, since the value for c\ can vary up to 1% when the 
lower limit is increased (the variation of C2 is smaller). We find, however, that the choice of 0.7 as the lower limit 
minimizes the reduced % 2 value of the linear fit. 
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k,=0 data, linear fit 



as 
-a 




FIG. 3. The derivatives of Sp with respect to the amplitude of the zero mode, evaluated at several values of 
a. The slope of this data is used to determine r max of the bilinear kernel Q(x — y), as explained in the text. 



make an ansatz for the kernel which imposes the finite range restriction on Q in a simple way: 



' yl, )xy '* yl " VmaX . (22) 

\x — y\> r n 



' max 



Given r max , Q{ki) is obtained by a Fourier transform of Q(x — y). To determine r max , we do a 
linear fit of the k^ = data 

1 fdSp 



(23) 

V \da J ao = a 

as shown in Fig. 3. Let the slope of the line be D. Then, from (21), 

c l -4c 2 Q(0))=D, (24) 
and we choose r max to satisfy this condition as closely as possible. We then have Q(ki) at all k. 

In Fig. 4 we plot the data shown in Fig. 1 together with the values computed for 

a{ l -c l -2c 2 Q{k L )) (25) 

(cf. eq. (21)) at a = 0.05. Agreement seems to be quite good in the entire range of k^. 

We have repeated this analysis at smaller volumes of spatial extension L = 12, 16,20. The 
results for c\ , C2, r max are shown in Table I. In Table II we record non-zero components of Q{x) = 
(—Vj)xO U P to |jc| < 3.2 and lattice volume 24 3 . For r max = 3, the last entry in the table should 
be replaced by 2(3, 1,0) =0. The rest of the non-zero elements of Q are obtained from the table 
via permutation symmetry, Xj <H- xj, and reflection symmetry Xj — > —xi, among the coordinate 
components. 
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L=24 data and finite range kernel 
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FIG. 4. A test of eq. (21) at a = 0.05. The derivative data of Fig. 1 is plotted against the conjectured fitting 
function a{\c\ — 2c2<2(&l)) with r max = 3. 
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Cl 


I"max 


12 


4.364(6) 


0.491(1) 


3.2 


16 


4.417(4) 


0.498(1) 


3.0 


20 


4.416(7) 


0.493(1) 


3.0 


24 


4.414(8) 


0.493(1) 


3.0 



TABLE I. Constants defining the effective Polyakov line action for pure YM theory, L 3 x 4 lattice, /3 = 2.2. 



X\ 


XI 


*3 


Q(x) 











2.38760 


1 








-0.22001 


1 


1 





-0.02357 


1 


1 


1 


-0.00774 


2 








-0.01279 


2 


1 





-0.00455 
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1 
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-0.00246 


2 


2 





-0.00160 


2 


2 


1 


-0.00111 


3 








-0.00200 


3 


1 





-0.00121 



TABLE II. Non-zero elements of the bilinear kernel Q(x) at r max = 3.2 and L = 24. 



B. Comparing the PLA to the underlying lattice gauge theory 

We now have a concrete proposal for the effective Polyakov line actions at various spatial 
volumes L 3 ranging from 12 3 to 24 3 , and which correspond to an underlying lattice SU(2) gauge 
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Polyakov line correlators, L=12 



Polyakov line correlators, L=1 6 
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(c) L=20 
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(d) L=24 



FIG. 5. A comparison of the Polyakov line correlation functions G(\x — y\) = {P x Py} as computed via 
lattice Monte Carlo simulation of the underlying gauge theory on a L 3 x 4 lattice at coupling j3 = 2.2, and 
via Monte Carlo simulation of the corresponding effective action Sp of eq. (19). Lattices are of spatial 
extension L = 12, 16,20,24 lattice spacings. Note that off-axis displacements are included. 



theory at /3 = 2.2 on an I? x 4 lattice volume. The actions are specified by eqs. (19), (22), and the 
constants in Table I. Above the lattice volume 12 3 , these actions are about the same. 

The crucial question, of course, is whether these proposed Polyakov line actions are correct; 
they are certainly different from the action suggested in ref. [6], which was derived for gauge 
configurations in a rather unrepresentative region of configuration space. There is one obvious 
and essential test: Do the derived Polyakov line actions reproduce the Polyakov line correlator 
calculated in the corresponding lattice gauge theory? Thus we compute, via numerical simulation 
of the Polyakov line action at L = 12, 16, 20, 24, 

G(\x-y\) = {P x P y ), (26) 

and compare the result to the same observable obtained from standard lattice Monte Carlo at 
/3 = 2.2 on an L 3 x 4 volume. The result for these four cases is shown in Fig. 5. Note that off-axis 
displacements are included, with xyz-components of the displacement x — y in the range [—4,4]. 

The data in Fig. 5 is limited to displacements of magnitude R < 7 lattice spacings, and it is 
interesting to consider larger displacements for larger lattices. This calls for higher statistics. In 
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On-axis Polyakov line correlators, L=24 




R 



FIG. 6. A high-statistics comparison of the Polyakov line correlation function G(\x-y \ ) = {P x Py) computed 
for the lattice gauge and effective theories, for displacements x — y parallel to the x,y or z-axes, and spatial 
volume 24 3 . 

Fig. 6 we compare the results for G(x — y) obtained on a 24 3 lattice for the effective theory, and on 
a 24 3 x 4 lattice for the SU(2) gauge theory, again at /3 = 2.2. In this figure we show the results for 
displacements x — y parallel to any one of the coordinate axes. The Liischer-Weisz noise reduction 
method [13] was used in obtaining the Polyakov line correlator in the lattice gauge theory, while for 
the effective Polyakov line action the correlator was obtained from 38,400 configurations (about 
two orders of magnitude more than was used in Fig. 5). 

The agreement of the correlators in the PLA and the underlying lattice gauge theory seen in 
Fig. 6 is extraordinary, and it persists down to magnitudes of order 10~ 5 . 2 While this is not a 
proof that Sp is the correct effective action, it is difficult to believe that agreement of Polyakov line 
correlators to this level of precision is coincidental. 

IV. POLYAKOV LINE ACTION FOR AN SU(2) GAUGE-HIGGS SYSTEM 

We now add a scalar matter field in the fundamental representation of the gauge group, thereby 
breaking explicitly Z2 center symmetry. The simplest case is a fixed-modulus Higgs field, and for 
SU(2) gauge theory the action can be written in the following way: 

S = P £ ^TrlUUUW] + K^Trl^(x)U^x)$(x + jl)} , (27) 

plaq L x,n Z 



2 The Polyakov line correlator derived from the Inverse Monte Carlo method in ref. [8] was displayed on a linear, 
rather than logarithmic, scale, and hence the precision of agreement with lattice gauge theory, in that approach, is 
difficult to judge. 
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where 0(x) is SU(2) group- valued. The work of Fradkin and Shenker [14], itself based on a 
theorem by Osterwalder and Seiler [15], demonstrated that the Higgs region and the "confinement- 
like" regions of the /3 — fc phase diagram are continuously connected. Subsequent Monte Carlo 
studies found that there is only a single phase at zero temperature (there might have been a separate 
Coulomb phase), although there is a line of first-order transitions between the confinement-like and 
Higgs regions, which eventually turns into a line of sharp crossover around /3 = 2.775, K = 0.705, 
cf. [16] and references therein. At /3 = 2.2 the crossover occurs at K ~ 0.84, as seen in the plaquette 
energy data shown in Fig. 7. There is also a steep rise in the Polyakov line expectation value as K 
increases past this point. 



SU(2) Gauge-Higgs, P=2.2, 16 4 Lattice 
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FIG. 7. Plaquette energy vs. gauge-Higgs coupling K at fixed j8 = 2.2, for the SU(2) gauge-Higgs theory 
with fixed Higgs modulus on a 16 4 lattice volume, showing a sharp crossover at K 0.84. 



We will work at /3 = 2.2 on a 24 3 x 4 lattice volume, but this time at Higgs coupling y = 0.75, 
which places us in the "confinement-like" phase a little below the crossover point. For these 
parameters, the Polyakov line has a VEV of (P x ) = 0.0515. Once again, we generate sets of 
fhermalized Polyakov line holonomies, and compute L^dSp/da^ as explained in the previous 
section. 

The derivatives dSp/da^ at = a are computed as before, and at each k^> the results are 
simply proportional to a. The constants c\,C2 are again extracted by a linear fit to the ki > 0.7 
data. However, the data at ki = is not strictly proportional to oc; there is also an a-independent 
constant contribution to the data. This fact can be seen in Fig. 8. The straight line is a best fit to 
LT^dSp/dciQ evaluated at ao = a for a in the range [—0.2,0.2]. The y-intercept of this line does 
not pass through zero, but rather through y = 0.0236( 14) . This implies that Sp must contain a term 
which is linear in P x , i.e. 

Sp = c J> x + X -c x £/^ 2 - 2c 2 Y,P*Q(x ~ y)Py , (28) 

x x xy 

and it is clear from inspection that cq must equal to the y-intercept in Fig. 8. It is also clear that 
only the k^ = mode contributes to the linear term, and is therefore invisible in the derivatives of 
Sp at ki > 0. We define Q(x — y) again by (22), with r max determined as in the pure-gauge theory. 
The final set of parameters for the effective Polyakov line action is given in Table III, and we plot 
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gauge-Higgs k=0 data, L=24 



gauge-Higgs k=0 data, L=24 
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linear fit 









-0.04 



-0.02 





a 



0.02 



0.04 



(b) close-up 



FIG. 8. The derivatives of Sp with respect to the amplitude of the zero mode in the gauge-Higgs theory, 
evaluated at positive and negative values of ao = a. (a) shows the full range of the data; (b) is a closeup near 
a = 0. The v-intercept of this data is non-zero, and determines the coefficient cq of the linear, Z2-symmetry 
breaking term in the effective PLA (28). 



the ki > data, together with the quantity 

a{^ Cl -2c 2 Q{k L )) vs. k L (29) 

in Fig. 9. 



L 


CO 


c\ 


C2 


r max 


24 


.0236(14) 


4.441(9) 


0.501(1) 


3.2 



TABLE III. Constants defining the effective Polyakov line action for gauge-Higgs theory, 24 x 4 lattice, 
J8 = 2.2, fc = 0.75. 



As in the pure gauge theory, the crucial test is to see whether the Polyakov line correlator (26) 
found from numerical simulation of the gauge-Higgs theory (27) agrees with the same observable 
computed in the derived Polyakov line action (28). The results are shown in Fig. 10. In this case 
the agreement between the lattice gauge Higgs correlator (black diamonds) and the correlator of 
the effective action (blue triangles), while fairly close, is not perfect. However, the result for the 
effective action depends very sensitively on the value of cq, and of course there is an errorbar 
associated with this quantity. In our best fits, co = 0.0236(14). With a little trial and error, one 
can find a value of cq for the effective effective action such that the corresponding correlator (red 
circles) agrees almost exactly with the gauge-Higgs value. This happens at a value cq = 0.02165 
which is not far outside our errorbars, about 1.4 a away from cq = 0.0236. 



V. CONCLUSIONS 

Motivated by the well-known sign problem, we have applied the relative weights method to 
determine the effective Polyakov line action Sp for both pure and gauge-Higgs lattice SU(2) gauge 
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L=24 gauge-higgs data and finite range kernel 
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FIG. 9. Same as Fig. 4 for the gauge-Higgs theory. We plot the data for the derivative L dSp/da^ vs. Rl 
against the conjectured fitting function a{^c\ — 2c2<2(^l)) with r max = 3.2. 



On-axis Polyaov line correlators, L=24 




I i I i I i I i I i I i I I 

4 8 12 16 20 24 
R 

FIG. 10. A comparison of the Polyakov line correlation functions G(\x — y\) = {P x Py) as computed via 
lattice Monte Carlo simulation of the underlying gauge-Higgs theory (black diamonds) on a 24 3 x 4 lattice, 
at couplings j3 = 2.2, K = 0.75, and via Monte Carlo simulation of the corresponding effective action Sp 
of eq. 28 (blue triangles, cq = 0.0236). Also shown is a simulation of the effective action with a slightly 
different value of cq = .02165 (red circles). 



theory. This effective action turns out to be a remarkably simple expression, which is bilinear in 
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the Polyakov line variables P x 



S P 



c J>, + \c, Y^Pl - 2c 2 Y,P x Q(x - y)P y 




Q(x-y) 



(30) 



\x-y\>r max 



with co = in the pure-gauge theory, and non-zero in the gauge-Higgs theory. Our results so far 
have been obtained at lattice coupling /3 = 2.2, and N t = 4 lattice spacings in the time direction. 
The effective action has been checked by computing Polyakov line correlators in both the effective 
theory and the underlying gauge theory, and we have found that these correlators agree quite well 
with each other. This is especially true in the pure gauge theory, where agreement persists down to 
correlator values of order 10~ 5 . These results, together with previous checks in the case of strong 
coupling [6], inspire some confidence that the method works. The next immediate step will be 
to understand how the couplings of the effective theory evolve as a function of coupling /3 and 
temperature l/N t . 

In the longer term our interest is the sign problem, and to address that problem it will be 
necessary to derive the effective Polyakov line action corresponding to lattice SU(3) gauge fields 
coupled to matter at zero chemical potential. In lattice gauge theory fixed to temporal gauge, the 
chemical potential /i is introduced via the replacement 



It is not hard to see that, to all orders in the strong coupling + hopping parameter expansion, the 
effective PLA obtained at ji ^ is related to the action at /x = by a simple substitution 



We will assume that this identity holds in general. The strategy is then to apply one or more 
of the methods [1-4], which were developed for solving Polyakov line actions with a chemical 
potential, to our derived effective action. If the sign problem is tractable in the effective Polyakov 
line theory, as suggested by the earlier work cited above, then it may be possible to extract useful 
results regarding the QCD phase diagram. 
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Appendix: Dynamical Fermions 

In this appendix we will just sketch how the relative weights algorithm can be applied in the 
case where there are two mass-degenerate dynamical fermions coupled to the gauge field. 



U (x,t = 0) -+e 



''^Uo(x,t = 0) , U%(x,t = 0) -»■ e~ Nttl Ul{x,t = 0) . 



(3D 



S^°[U X , U*] = S^ =0 [U X -)• e N >»U x , Ul -)• e^U*] . 



(32) 
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Taking D as the Dirac operator, and M = D^D, the integration measure is 

e Sl = detM e Sw , 

where Sw is the Wilson (or other improved pure-gauge) action. Defining 



M 



(m) 



M 



U (x,0) = ut } ,U k (x,t 



(A.l) 



(A.2) 



we have 



exp[As£" +1) 



7 

I detM( m+1 ) 
\ detAfW 



exp[A5'; ,+1 



(A.3) 



where (...) m signifies, as before, the expectation value in a probability measure exp[5^]/Z m . 
The VEV can be evaluated via hybrid Monte Carlo, but in this case there is the question of how 
to evaluate the ratio of determinants in (A.3). In the most straightforward approach, defining 

5M ( m +l) = M {m+l) _ M (m) and writing 



detM ('"+i) = det|M(' n )(l + [mW]- 1 ^ 1 ))} 



(A.4) 



it is not hard to see that 



detM('" +1 ) 
detMH 



exp pTrlogjl + [M^]~ l 8M^ m+V ' 



exp[Tr{[MW]- 1 5M( m+1 ^ 



Therefore 



exp [AS 



'm+D 



exp 



Tr\[M^]- l 8M {m+V ' 



e w 



(A.5) 



(A.6) 



It is best not to evaluate the trace Tr{[M( m )] l dM^ m+ ^} directly, because of the computational 
expense, but rather indirectly, by the method of "noisy pseudofermions" (c.f. section 8.4 in ref. 
[17]). 

A way of completely avoiding M _1 in the observable is to introduce another set of pseud- 
ofermions (p, distinct from the pseudofermions used by hybrid Monte Carlo, so that we may write 



detM(" 7+1 ) fDq>exp[-<p*MW(p] 



Defining 



we have 



detMH /D<pexp[-<ptM(' n+1 )<p] ' 

M (m+1) = ( pt M (m+l) (p _ (p t M (m) (p ^ 

detM( m+1 ) 



detAf( fi 



exp[AA(" 1+1 )] 



m+l 



(A.7) 



(A.8) 



(A.9) 
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where ((...)) OT+ i refers to the expectation value in the measure proportional to 

expt-^Ar^V] . (A. 10) 

The final result for the relative weight is 

exp[A5^ +1) ] = ^ ((exp[M( m+1 )])) TO+i exp[A5i; i+1) ]^ . (A.ll) 

Again, the expectation value (...),„ would be computed via the usual hybrid Monte Carlo algo- 
rithm. Eq. (A.ll) is to be used, as before, to compute path derivatives, hopefully leading to an 
expression for the effective PLA at \x = 0. The PLA with a finite chemical potential would then 
be obtained from the substitution (32). 

Apart from computation cost, we believe that the addition of dynamical fermions does not pose 
any problems of principle to the derivation of the effective Polyakov line action via the relative 
weights approach. 
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